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Abstract 


With numerous foundry services being made available for implementing 
systems on VLSI chips and with the rapid stndes that are bemg made m the 
fabrication technology of these chips it is imperative from the point of view of 
economics to be able to migrate any design from one foundry specific technology to 
another as well as from a present generation fabrication technology to the next 
generation fabncation teclinology within the same foundry As the ability to integrate 
and pack more devices within a die of silicon increases with every new generation of 
the fabrication technology the complexm of digital systems realizable in a single 
chip has also grown by leaps and bounds To circumvent implementation bottlenecks 
a lot of research work has been carried out to fully automate the design of VLSI chips 
under prespecified area and performance cnteria The focus of the present thesis is on 
an important aspect of VLSI chip design known as Physical Layout Design 

The basic issue addressed here is the following Is it possible to port and 
reuse existing cells in a particular generation of a fabncation technology to a new, 
but evolving generation by using the mask layout descriptions of the existing 
cells As the creation and validation of the mask layout descnption of cells m anv 
fabrication technology is a time consuimng error prone tedious and costly process it 
IS important to be able to make the best possible use of the existmg layout resources 
accumulated from the imtial fabrication process This process is known as technology 
migration 
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In this thesis study of some existing software tools and algonthms that ha\e 
been employed for technology migration has been earned out Some new approaches 
and the correspondmg algonthms for carrying out technology migration has been 
proposed These algonthms ha\e been incorporated in a technology migrator called 
TECHMIG TECHMIG has been de\ eloped as a part of this thesis TECHMIG 
takes as its input the mask lav out descnption of an existing cell expressed m the 
industry standard Caltech Intermediate Format (CIF), along with a set of design 
rules specific to a new fabncation technology or a new foundry and a set of user 
specific design constraints It then retargets the imtial cell by producing its mask 
lav out description in CIF which is appropnate for the targetted fabncation technologv 
or the foundry 
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Chapter 1 
Introduction 


State of art VLSI fabncation technologies allows fabrication of devices with 
geometry smaller than 0 25 micron Advances in VLSI fabncation will further lower 
this \alue in the near future With smaller geometry the ability to fabneate systems on 
chip IS already being demonstrated by many design houses With the possibility of 
integrating increasingly complex designs in a single die of Silicon chip it is 
imperative to improve both design productivity as well as design quality by 
leveraging existing standard cell libranes developed for commercial purposes and 
research work A large component of any library developmental work involves 
physical design and layout of cells or circmt elements for a presenbed set of foundry 
specific design rules The library development work is in general earned out by a very 
skilled set of layout designer With typical libranes that are bemg made available by 
foundries having more than a few hundred cells and megacells any library 
development for new VLSI fabncation processes and technologies with newer design 
rules IS by and large a very time consummg error prone and tedious activity Each 
new library development generates very few newer cells but retargets existmg cells m 
older technologies for the new design rules Since the design rules are foundry 
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specific a cell designed for a particular foundry is useless for another foundry which 
obeys an entirely different set of design rules If we have to fabncate the same cell in 
the second foundry we have to completel> redesign the cell 

Generall> renewal of VLSI fabncahon process is an evolutionary process 
with incremental changes bemg mtroduced in an existmg set of design rules As such 
a set of layout resources accumulated from an old fabncation process must be best 
utilized to gain full advantages of the newer fabncation process In this thesis we 
address this issue Specificallv we address the following question js it possible to 
use the full masked version of a cell designed in an older technology to create 
another cell with identical functionality and resource usage, while maintaining the 
saint layout topology that was present in the original celP This process of exploiting 
the reusability of layout resources and translating it from one set of design mles to the 
other IS called "Technology Migration" 

fo improve the design quality and design productivity every manual process 
must be replaced by a well organized set of automatic/semiautomatic tools Most 
design houses have a history of physical libranes that have been developed for 
different foundries and these serve as the starting point of the retargettmg work These 
cell layouts are manually edited and migrated to the target design rules set This task 
IS often very time consuming and manpower intensive The motivation for automation 
arises due to the following factors 

• As fabncation technology advances rapidly it often results in a changed 
fabrication process The new changes result in a large scale change in the 
design mles 

• Circuit performance based on submicron technology is greatly affected by the 
signal delays Hence there anses a technical issue as to how to optimize shapes 
and sizes for different kmds of layout elements 
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« In order to reduce the turn around tune for the custom VLSI design it is of 
practical importance to be able to provide m advance a full set of cell libraries 

The present thesis is devoted to addressmg this problem and to de\elop a 
highl\ automated CAD tool for migraung cells from one foundry to the other In this 
thesis we restrict ourselves to the problem of technology migration for designs based 
on the standard cell approach In this chapter we first descnbe in details the process 
of teclinology migration and the basic rmgrator model that we have developed 


1 1 Migration Tools, An Overview 


There are two broad categones of methodologies capable of efficientlv 
letargetting layout between vanous processes These are 

• To generate layout inside tools that can retarget the results based on different 
design rules 

• To migrate GDSII or CIF libraries from one process to another 

Each of these methodologies have their advantages and shortcomings In both 
the methodologies once the migration technology file is written for one process the 
second process migration will be verv fast assuming an evolution and not a 
revolution in the new VLSI fabncation process Layout migration is inappropnate 
when there are drastic differences between the source and the target process e g 
migrating from a 2 metal layer process to a 3 metal layer process Following are some 
of the characteristics of the two different approaches of technology migration cited 
above 
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Generating a new layout Here e\er> cell should be generated inside the tool The 
user needs to be an expert m not onl\ in laving out the cells but also in understanding 
the specific tool philosophy The final results are not perfect in either methodolog\ 
but m this case the user has to work within the tool to get the desued results The 
dement of this approach is that once out of system back annotation for subsequent 
improvements to the source is not possible Further it can not be used to migrate old 
libraries or chips that were not previouslv de\ eloped with this methodologv The 
principal advantage with this approach is that cells/blocks topology can be altered 
inside the tool 

Migrating the GDSII or GIF In this approach one can generate the libranes with 
im kind of layout editor or else one can use old cell layouts These tools use the 
standaid GIF or GDSII formats as input and the user may not ha\e to be a la\out 
expert to use the conxerter If the result of the conversion is not satisfactory a normal 
la\ out editoi can be used to change the source The principal shortcoming of this 
approach is that the cell topology cannot be changed to take the advantage of a 
different piocess We adopt this approach tor our technology migration tool 

Migration tools also fall into two other categories viz conversion and 
compactor tools A conversion tools converts a layout from one process to another A 
tool such as Dracula from Gadence Design Systems can do the conversion, but 
needs additional programming to do compaction 

By contrast a compaction tool accepts an existing layout in form of GDSII or 
GIF format then applies the design rule files for a new process The compactor then 
shrinks the existing layout to a smaller geometry layout converting the layout m the 
process In operation the compactor shrinks the cell in honzontal direction until a 
design rule is broken then performs a similar task in the vertical direction Afterwards 
It begins to apply additional algonthms to further reduce the cell size We follow the 
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compaction approach rather the conversion approach for migratmg the layout In this 
section we describe the basic technology migrator tool TECHMIG developed as a 
part of this thesis work its user mterfaces and its internal structure In subsequent 
chapters we will descnbe the process of carrying out layout compaction 


1 2 TECHMIG 

12 1 User Interfaces 


The various user intertaces of the TECHMIG tool is shown in Figure 1 The 
piincipal component is the central block TECHMIG The other blocks are the user 
intertaces These are described as follows 

12 11 The Input Layout 

The b isic input <-o the technology migrator is a descnption of the ongmai cell to 
be migrated Almost ail the compaction algorithms reported in literature take a 
svmbohc description of the layout However in order to capture the full mask features 
ot a layout the symbolic descnptions are mconvement Poor to fabncation the cell 
must be available in the mask layer descnption format rather the symbolic format The 
svmbohc representation is mconvement and infeasible m representing some extended 
lemons such as h well and p-well in CMOS technology Thus shape descnption at 
mask level rather than the symbolic descnption was the preferred choice as an mput 
to the tool We assume the standard Caltech Intermediate Form (CIF) representanon 
of a layout The details of CIF format are given in [1] However we deviate shghtly 
to tollow the Compass tool compatible CIF version 



figure 1 The technolog^ migrator TECHMIG with its user interfaces 


12 12 The Design Rules 

Layout rules or design rules are constraints imposed on the geometry of an 
integrated circuit layout in order to guarantee that the circuit can be fabricated with an 
acceptable yield These rules are essential for any design implementation to match the 
fabrication process specifications Adhering to the design rules is essential for anv 
product to be manufacturable and reliable after fabrication The mam objective of a 
set of design rules is to ensure optimum performance with minimum implementation 
area without compromising production yield reliability or functionality of the system 
being integrated into a chip The design rules act as an mterface between the circuit 


6 







designer and the process engineer Circuit designers m general want a tighter smaller 
design that results in higher performance and higher circuit densitv The process 
engineer on the other hand wants a reproducible and high yield process In general 
design rules represent the best possible compromise between mtnnsic dewce 
performance and vield In general more conservabve the rules are the more Idcelv the 
circuit will function but with relaxed performance On the other hand if the la-sout 
rules are verv aggressi\e it results m better circuit performance due to reduced 
parasitic and improved charactenstics The design rules primanly address two issues 

• The geometncal reproduction of features that can be produced bv mask 
making md the lithographic processes 

• The interference constraints between different layers 

The tundamental uniU in the defimtion of a set of design rules is the 
minimum line niclt/i which is defined as the minimum mask dimension that can be 
said) iranstened to the semiconductor matenal It is set by the resolution of the 
putcinmg piocLSS At least two basic approaches are taken m descnbmg the 
dimensions in design rules 

Scalcable Dimensions The scaleable design rules also called the lambda (X) based 
rule was made popular by Mead and Conway [1] It defines the entire set of design 
rules as a function of a single parameter called lambda (X) which is specific to the 
particular tabucation process In Mead and Conway [1] the basic umt of length 
measurement is equal to the fundamental resolution of the process itself Lambda is 
defined as the distance by which a geometncal feature on any one laver may stray 
from another geometncal feature on the same layer or on another laver considering 
all the processing factors and adding an appropnate safety factor The value of 
lambda is set by phenomena associated with the fabncation process such as over 
etching misalignment between mask levels distortion of the silicon wafer ( run out ) 
due to high temperature processmg and overexposure and underexposure of resist 
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In this method the rules are chosen so that a design is easily ported over a cross 
section of industnal process Scalmg of the mimmum dimension is accomplished bv 
simpl\ changing the value of lambda For a given process 1 is set to a specific 
absolute value and all design functions are consequently translated into absolute 
numbers The minimum line width of the process is set to 2X. For instance for 1 2 
micion process 2 equals to 0 6 micron The detail descnption of lambda based rules 
is j;,iven m [1] and [2] 

Absolute Dimuisions (microns) Here all the design rules are described as its 
ibsolute \.alue For example the minimum width of active diffusion area mav be 
specified as 1 micion Almost all the industries follow this stvle than following the 
1 unbda based lule The reasons behind the adaptability of micron based rule are 

• Adopting the lambda based rule exploits the advantages of linear scaling 
of the ciicuit elements However as the design complexm increases with 
simultaneous shrinkage of device size linear scaling can lead to design 
rule violations Linear scaling is only possible over a limited range of 
dimensions (for instance between 3 micron and 1 micron) When scaling 
o\ei larger ranges (for instance into sub micron range) the relations 
between different lavers tend to vary in a nonlinear wav that cannot be 
adequately covered by linear scaling rules 

• Scaleable design rules are conservative As they represent a cross section 
over different technologies they have to represent the worst case rules for 
the whole set which results m an over dimensioned and less dense design 

For these reasons scaleable design mles are avoided by mdustnes As circuit 
density is the pnmary goal in industrial design most semiconductor mdustnes tend to 
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use micron rules Ho\\e\er the scaling and portmg of designs between technologies 
under these rules is more demanding and has to be done either manually or bj usmg 
rule The reasons behind the adaptability of micron based rule are 

T\PES OF DESIGN RULES 

At a sutficienth high conceptual le\ el all CMOS processes use the following features 

• Two diffeient substrates and or wells (CNW CPW) 

• Active regions tor both p and n transistor as well as substrate taps (CND CPD) 

• One or more pol\ sillicon la\ ers which serves the purpose of gate of the transistors 
IS well as the interconnection layer (CP) 

• One or moie levels ot metal interconnection layer fCM CM2) 

• Contacts and \ las to provide mterlaver connections (CC CC2) 

Detailed description of these lavers axe given in [3] The design rules reflects 
the interaction ot these 'avers in the lav out plane Broadly we can classify the whole 
set oi design lules into following categones 

• Size Rules 

• Separation Rules 

• Overlap / Extension Rules 

Si^.e Rules These rules refer to the minunum width of a wire on a particular laver 
This IS imposed to prevent thin features such as wires being narrow enough to break 
the connection and also to prevent electromigration The minimum feature size of a 
device or an interconnect is determmed by the Ime patterning capabihty of 
lithographic equipment used for IC fabncation Interconnect Imes usually run over a 
rough surface unlike the smooth surface over which active devices are patterned 
Consequently the minimum feature size used for mterconnects is somewhat larger 
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than that used for actne de\ices The vanous design rules under this categorv are 
gi\en in [3] Figure 2 gives an illustration In Figure 2(a) the rmnimnm required 
width of a diffusion laver is shown and m Figures 2(b) and (c) imnimum width for 
poly and metall la\ er are shown 

Separation Rules Different features on the same layers or on interacting layers need 
to be separated from each other These separation type design rules prevent the 
separate and isolated features from accidentally short circuiting with each other 
Figure 2 illustrates some of the spacing rules These spacing rules can also be 
specified conditionalK For example the spacing between two metal layers depends 
upon the thickness of the wires and the extent to which the two wires run in parallel 
The lequired spacing between two wide metal wires is somewhat larger than the 
rcquiied spacing between two narrow metal wires Similarly the required spacing 
between two parallel and long wires is somewhat higher than the required spacing 
between two short wires These types of rules which depend upon certain 
chaiacteiistics of the primitives are called the conditional rules TECHMIG has the 
capability of supporting such conditional rules A detailed listing of various spacing 
lilies both conditional md unconditional is given in [3] The separation rules for 
three different lavers are shown m Figure 2 



(a) Diffusion 


(b) Poly 


(c) Metall 


Figure 2 Illustrating size and separation rules for diffusion poly and metall layer 
All dimensions are m hundreds of microns 
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Overlap/Extension Rules The overlap rules ensure the integrity of the topoloeical 
features m a lavout Tlie extension rules particularly refers to the extension of polv 
region bevond tlie diffusion region delSning the transistor This pre^ents the 
formation of incorrect transistor or short circuited channel In Figure 3(a) a contact 
overlap rule is shown and in Figure 3(b) the extension rule for a transistor is shown 
Detailed listings of the vanous o\ erlap/extension rules are given m [3] 


CM 

60t 

*60 


CP 


(a) 





^60^ 


60 




CND 



(b) 


1 igure j Illustrating overlap and extension rules (a) CC CM and CC ~ CP overlap 
iLile and (b) CP CND extension rule Shaded area is the transistor All dimensions 
ire in hundreds of micron 


REPRESENTATION OF THE DESIGN RULES 

Specilic to a fabrication process the number of design rules typicallv vanes 
fiom one to a few hundred Tliese are normally specified in plain English bv a one 
sentence description Additionallv graphic illustration may be provided for the vanous 
cases ot applying these rules The common representation of fabncation rules is by a 
matnx of simple spacing rules between pans of masks with size and overlap rules 
stored in separate tables Presently vanous design houses follow the EDIF format for 
design rule descnption However we follow a uniform framework of design rule 
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description as descnbed in [4] Here we bnefly descnbe the format for design rule 
description The followings are some of the chaxactenstics of this method of 
representation 

• All rules are formulated m the same maimer mcludmg space rules size 
rules o\ erlap rules and extension rules This is essential for the compaction 
of shape objects 

• Conditional rules are allowed We support the conditional rules related to 
some specific topology Provisions are not made to represent other 
conditional rules related to parasitic resistance capacitance and reliabilitv 

• A. linked list structure instead of a full matrix form is used for the lav out 
lilies to i educe the space complexity In matrix form we also have to reserv e 
spaces tor nomnteracting hvers for which no design rule exjst 

The boundaix of a mask shape is a set of lines with each line having two 
sides namely imidc or outside The inside refers the side of a line towards the mask 
ind the outside is the side not included in the mask Since the rules depend upon both 
the mask level and the side of the boundarv an edge is defined as a boundarx line 
with a side attribute That means 

edge == T(mask side) 

wheie T is the transformation applied to the mask depending upon the side attnbute 

The required transformation T is defined as 
T(mask Inside) = mask 
T(mask Outside) = -mask. 

For example on an n diffusion laver (CTW), CND represents an mside edge while 
CND represents an outside edge 
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A simple rule between a pair of masks can be expressed as edgel edge2 
Value where Value is the absolute value of the design rule between the 
corresponding edges in hundreds ot micron (standard umt of the dimensions m GIF 
representation of the la\out) Some examples of simple rules are shown m Figure 4 In 
Fij.ure 4(a) the overlap and size rules are specified and m Figure 4(b) the spacing 
rules are specified Here maskl is a contact layer and mask! is a metal 1 layer 

Tlie size mle for the metal wire is given as CM CM Wmin Wmm bemg the 
minimum required width of the metal wire If the minimum overlap of a metal layer 
on the contact is Min_Overlap the corresponding rule is specified bv CM CC 
\Iin_Overlap Similarh the minimum spacing rule between two poly la\ers is 
specified as CP CP Min_space 

Conditional rules are represented by augmenting the edge descnption with 
tnothei attribute For example the lule between two wide metal wires on the metal 1 
liver can be represented bv -CM_W CM_W value_w and the rule between two 
11 urow metals similarly can be represented by CM_N CM_N value_n 

flic detailed listing ot vanous upes ot design mles under these categones (size 
spice overlap) IS specified in [3] 


\ 

CM 

► 

CC 


(b) 

Figure 4 Illustrating the representation of design mles 


maskl 


mask2 


(a) 
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12 13 User Constraints 


These refer to some special types of rules as per the requirement of the user 
commonl} known as the user constraints (UC) These include the power bus 

width connector positions transistor sizes cell height etc pertinent to the 
cell The fourth t\pe of rules le the user constraints does not depend on the 
particular fabrication process rather are specified by the user dependmg upon his/her 
requirements These being \er> specific to a lavout are represented in a separate 
format 


The mles regarding transistoi resizing are stored in a table called Transistor 
Rcsuing Table (TRT) which has information regarding the position of the transistor 
m the original la\out it s fV/L \aiue m old technology and its W/L \alue m the new 
technology The other user constraints regarding power bus (VDD and USSj width 
cell height and connector placement are provided bv the user at run time 

12 14 Output Lavout 

As the tool does not change the topology of the ongmal layout the 
output layout has the same format as the mput layout Thus our tool produces the 
output cell m CIF format compatible with Compass CIF version 
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1 2 2 Internal Structure of TECHMIG 


The different blocks that tire tool consists of are shown m Figure 5 These are 
described as follows 



1 DR Design Rules 2 UC User Constraints 

j CiF_IN Input GIF description ot cell 4 I P Input Processor 
0 0 P Output Processor 6 C G G Constraint Graph Generator 

7 C G S Constraint Graph Solver 8 L O Layout Optimizer 

Figure 5 Illustrating the internal structure of the technology migrator TECHMIG 


15 





12 2 1 Input Processor 


Tlie input to this block is a Compass tool compatible CIF description of the 
source layout the design rules files and vanous user constramts The mput processor 
(I P) reads these information processes it and stores the relevant information in some 
internal databases like an arrav of polvgons or a linked list etc The representation of 
design rules has alreadv been described and these are stored in a linked list format 
Here we give a brief descnption of the internal representation of th.e input cell The 
/ P handles this task of conv enmg the CIF description of the layout to this internal 
representation 

Li} out patterns tor functional cells are composed of layout elements Each 
element is lepresented as a pohgon with attributes such as shape size location and 
the mask laver (such as diffusion metal etc) Such an individual element is also 
icfeiTed to as a primitive We use the term polygon instead to refer to either a 
primitive or an clement 

In almost all published articles the basic layout pnmitive is restricted to only 
icctangles An^ mask polygon of an arbitrary shape in the Manhattan geomem is 
sliced into a number of rectangles In our approach we overcome this limitation and 
instead represent the polygon as a single entity Tins form of representation helps m 
manipulating structures with non Manhattan geometry (e g 45 degree bent wires) and 
results in a constraint graph with minimum number of nodes This accelerates the 
subsequent steps of involving constraint graph solution and optimization The second 
adv antage is due to the fact that every edge of the polygon contnbutes to a node m 
constraint graph so the number of nodes m constramt graphs due to a particular 
polygon is equal to the number of edges in the polygon However if we take the 
rectangles as the basic layout entit> and slice the polygon mto a number of rectangles 
the number of nodes conmbuted to the constramt graph depends upon the slicmg 
pattern as a particular polygon can vield many different combmation of rectangles for 
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different slicing patterns Optimal slicing of polygon is not well studied in literature 
and slicing to rectangles is not possible for layouts with non-Manhattan geometiy 
The entire layout is represented as a linked list of polygons The constraint generation 
step utilizes this linked list representation of the layout to analyze all the aspects of the 
lav out like separation overlap identification of transistors etc 

1 2 2 2 Layout Compactor 

A compactor takes as its input a VLSI lavout and produces as its output an 
equivalent layout ot smaller area The compactor moves components and wires in the 
hvout plane to optimize for two basic goals that the layout should be of as small area 
is possible and that it should be design rule correct The compactor moves subcells 
onlv in the plane preserving the cell topology Auxiharv functionalities such as 
vvirclength minimization automatic jog insertion contact optimization etc are 
generally incorporated to improve its capability to get a better layout according to 
specified performance criteria The compaction problem is stated as a mmimization 
pioblem that is find the minimum of the area function subjected to linear and 
nonlinear constiaints 

The compaction algonthms can be classified as one dimensional and t)vo 
dimensional compaction according to the movements of components during the 
compaction process In one dimensional compaction the mdividual components are 
moved only m a single direction either X or Y in a smgle pass For a complete 
compaction of a flattened layout two passes of ID compaction are necessary In 2D 
compaction a single pass can move a component in both directions 

Compactors can also be classified accordmg to the algonthms used to position 
the components The major techniques of the compaction are the constraint graph 
compaction (CG) virtual grid compaction (VG) and zone refining (ZR) techmque 
The descnptions of these approaches are given m [5], [6] and [7] We use the 
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constraint graph approach as it is the most popular and widely used method. We 
describe in details this approach in Chapters 2 and 3 The zone refining technique for 
1C layout compaction is given in [8] 

Several one dimensional and two dimensional compaction algonthms are 
available in literature of which a versatile one is presented in [9] which is based on 
the solution of longest path In [10] the compaction problem is translated mto a 
mixed integer linear programming of a special kind and a graph based branch and 
bound algorithm is used to solve the problem The literature presented in [1 1] uses an 
aLorithm almost parallel to that piesented m [9] except for the extra provision for 
considering the grid constraints li anv The problem of 2D compaction has been 
shown to be NP complete in [12] In [12] a branch and bound optimization (BBO) 
method is proposed tor 2D compaction the objective function being the minimum 
area of the bounding lectangle The problem of automatic jog insertion dunng the 
compaction process is not well studied in literature 

We restrict ourselves to alternating passes of ID compaction and formulate 
the compaction problem as a lineai programming problem We then solve this 
problem to minimize the bounding rectangle area The total wirelength serves as the 
auxiliary optimization function In our approach the entire compacnon problem is 
split into the constraint generation constraint solution and optimization phase In 
block diagram notations these steps are shown in Figure 5 In the first step the 
various design niles are imposed on the imtial lavout This requires the construction 
of the constraint graphs for both the X. and the Y direction The constramt generation 
IS the most important step of the compaction process as it requires a detailed analysis 
of the layout features using computational geometry After imposing the entire set of 
design rules the optimization problem is modeled as a Imear programming problem 
Using this formulation we search for a solution which leads to a layout with 
minimum area and vvhich satisfies all the design rules Wire length minimization is 
performed as a post optimization phase to additionally improve the layout to meet the 
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desired performance specifications Wire length rnmimization reduces the net 
parasitic resistance ot the la\out bv decreasing the wire length m highly resistive 
lavers and increasing it in comparati\ ely low resistive layers 

Various algonthms for generatmg the constraints are descnbed m Chapter 2 
In Chapter 3 we d scuss the \anous aspects of constraint solving and wire length 
optimization and the corresponding, algonthms 

1 2 2 3 Output Processor 

In the compactor we translate the layout features and design rules into a graph 
called constraint graph and tnen solve the compaction problem using graph theoretic 
ilgoiitlims The hnal stage in t<.chnologv migration consists of output processing In 
this phase the ^^laphs lepresentin^ the positions of the lav out pnmitives after the 
optimisation phase are translated into the CIF format Thus in this stage we obtain 
the output cell suitable tor a target toundrv 


Chapter 2 


Constraint Generation 

2 1 Intioduction 

One ot the basic task in constraint graph based compaction is to translate the 
design lules into a set of constraints and imposing those constraints to the lavout 
dements These con'itraints reflect the sizes of the primitives their relative positions 
in the lavout and some specific design rule requirements pertinent to the target 
fabrication process Keeping all these factors in view we can broadlv classify the 
whole set of constraints into following categones 

• Constraints specifying minimum required width of a particular layer 

• Constraints to represent minimum required spacing between two interacti\e 
lavers 

• Conshaints to represent o\ erlaps and extensions between primitives 

• Transistor resizing constraints 

• User constraints to fix the cell height power bus width, and connector 
positions in the layout 
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In this chapter we descnbe vanous algonthms to generate these constraints 
When the entire set of constraints is represented by a graph the resultant graph is 
known as a constraint graph Two graphs one to represent constraints m honzontal(X) 
direction called the honzontal constraint graph (Gpj or HCG) and the other 
representing set of constraints in vertical (Y) direction called the vertical constraint 
graph (Gy or VCG) needs to be generated We descnbe these constraint graphs bnefly 
as follows 

The layout is flattened and represented in a two dimensional (XY) plane Each 
edge of a primitive eg a rectangle is represented as a node in constraint graph A 
vertical edge contiibutes to a node in HCG and a honzontal edge contnbutes to a node 
in VCG Any inclined edge is repiesented as a node in both Gff and ui G^ Two nodes 
in 1 constraint graph are linked bv an arc if there exists any design rule which relates 
the conesponding edges of the primitives m the layout Thus the two graphs G^^and 
G[ encode the set of constraints G// encodes horizontal constraints while Gy encodes 
\ ertical constraints Each vertical or inclined edge in the layout is a node x^ in the 
jjiaph G^ and each honzontal or inclined edge is a nodey, in the graph G, The graphs 
uc duected acjclic graphs Each directed arc in the graph corresponds to an 
inequalitv Each iic (\, Xj, W,^ his a weight W,j associated with it representmg the 
constraint between the vertical edges corresponding to ^'he nodes x, and Xj respectively 
Implicitly the arc lepresents the inequality Xj-X^ >W,j 

Similarly each arc (x,^j) is associated with a cost parameter C,j that roughly 
indicates the resistance of associated layer This cost parameter serves as a secondary 
objective tor the compaction problem Exact formulation of cost fiinction and cost 
parameters are explained in Chapter 3 

Formally we define G„ = pC E„W„CJ iiSB. directed weighted graph X={xJ 
is the set of nodes in the graph one for each edge of layout and two special nodes 
(source) and (target or sink) corresponds to the leftmost edge m any cell in the 
layout and corresponds to the nghtmost edge 
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E { <x„ Xj> I x, Xj sX} IS the set of arcs in the graph, one for each mequalit> 

I <x„ Xj> e E„ W,j^ € R} is the set of weights specifying the separation 
between the nodes v, and Xj 

E\ ~ ] <;x, x^> e C j e R} \s the set of cost parameters related to la\er 

resistance 



Figuie 6 Illustrating a spacin^, constraint is the minimum required spacing 
between primitives -4 and B 



(a) (b) 


Figure 7 Illustrating the minimum width constramts Here is the mimmum 
required width of the primitive A 
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Some examples ot constraint graph are shown m Figure 6 and Figure 7 In 
Figure 6 two pnmitnes A and B are required to be spaced by a distance of at least 
Figure 6(b) shows thi_ constraint graph Here the cost value is set zero In 
r igure 7 for a primitu e A its minimum width should be umts and the layer cost 
is C units This is represented in graph domam as shown m Figure 7(b) 


2 2 Minimum Width Type Constraint Generation 


In this section we describe the algorithms to generate constraints that relate to 
the minimum requiicd width ol an entit\ This is dependent on the type of the layer 
issocnted with the entitv The mimmum required widths for each laver type 
corresponding to certain fabrication process are given in [3] In software 
implementation this minimum width tvpe of design rules are stored in a table whose 
indiMdual enti\ has three fields 

• La\cr Specities the charactenstic layer for the mask pol\ gon 

• bpccities the minimum required width for that layer 

• Cost Specitics the sheet resistance ot the layer 

In this section we will desenbe the algorithms to impose the mimmum width 
type constraints to the lavout and to add those constraints to the constraint graphs 
and Gp Following procedures tmd the mimmum width and the sheet resistance of a 

particular laver 

Procedure MinJVidth (Layer) 

Begin 

Curr_entry = first entry of the min_ width design rule table, 

While (All entries not searched) do 

Begin 


23 



If (Curr_entry layer = Layer) return (Curr_entiy n 

Curr_entry = next entr\ of the minjvidth design rule table 
End while 
End 

Procedure Min_}pidth (Layer) 

Begin 

Curr_Lntr\ - first entry of the min_width design rule table 

While ( All entries not searched) do 

Begin 

If (Cm rjentry -^la\ei — Layer) return (Curr_entry Cost) 
Curr_entr\ = next enti\ of the mmwidth design rule table 
End while 

End 


22 1 The Scanline Algol ithm 


To add the minimum width type constraints we take one primitive at a tune 
Depending upon its layer we determine the required minimum width as per the 
pioctdure \Iin_Wulth( ) Then we detect the edge pairs of the primitive which are to 
be related bv a wadth constraint Detecting all edge pairs involves ruimmg a scan Ime 
thiough the primitive m a direction perpendicular to the direction of constramts e g 
hoiizontal scanline mnning in vertical (Y) direction for adding honzontal width 
constramts and vertical scanline lunning m honzontal direction to add vertical width 
constramts Tlie following is an illustration of the scanlme techmque to detect all edge 
pairs which are to be constrained honzontally 
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A scanlme X AT is run in +K direction in the la>out plane An edge of the 
polygon IS inserted to the scanlme when its bottom end pomt is encountered and 
deleted from the scanlme when its top end pomt is encountered After inserting all 
ed^es that corresponds to a particular position (Y coordinate) of the scanlme adjacent 
edge pairs are checked to find out whether the edge pairs are exactly to be constramed 
It an\ arbitrar\ point intermediate to two adjacent points on the scanlme remains 
within the polygon then the correspondmg two edges represent the opposite edge 
pans of some portions ot the pol\ ^on and hence have to be constrained 



Figure 8 Illustrating the scanlme technique for generatmg mimmmn wndth 
type constraints 

In Figure 8 a primitive A is shown which is to be constramed through the size 
rules m honzontal direction A honzontal scanlme X’X is run from the bottom to the 
top of the polygon Its positions at distmct mstances are shown At position Yg, the 
edges e, and e, are inserted to the scan Ime It is found that these two edges enclose 
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certain portions of the pol>gon in honzontal direction which means that these two 
edges are actuah> to be constrained bv the width constraints and hence a width 
constraint is joined between edges Cj and Then the scanline is moved to the 
position 7/ where it encounters the bottom end points of the edges Cj, e and 
The edges C; e , and are inserted to the scanlme Adjacent edge pairs e_ and 
tj and Cfi are constiained as the% enclose some portions of the polygon Howe\er the 
adjacent edge pairs e and are actually not to be constramed We take this decision 
b\ predicting a point intermediate to e and and checking whether that point 
lemains within the pohgon It it lies within the polygon then the edge pairs are 
ictuallv to be constiained b\ width constraint At position Y the bottom end points of 
edges e tj and c, are encountered bv the scanlme hence e, and are deleted 
horn the scanlme biniilarh at position Fj the edges ej are deleted from the 
scanlme Figuie 9 shows the constraint graph generated bv this algonthm A formal 
description ot oui ilgorithm to j^enerate size constraints in both honzontal and \ enical 
direction corresponding to a layout entit) is giaen below 


W 

mm 



Figure 9 The constraint graph showing Constraints only as genertated 

by scanlme technique 
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Procedure Add_ Min_ WidtIi_Constraints_to (entity) 

Begin 

^^nui ~ yfinimiun width of layer corresponding to entity, 
if (entity -^t}pe = box) 

Begin 

Find opposite edge pairs of the box Join the constraints of weight 
between the corresponding nodes in the constraint graph 
return 
End if 

if (entiU, ->t\pe= polygon) 

Begin 

Step 1 Cteate an arra\ Hf] whose elements are the end points of the 
horizontal and inclined lines Similarly create an array VfJ whose elements 
are the end points of vertical and inclined lines 

Step! Sort the elements of H[J array according to the following strategies 

• Take Y coordinates of the points as primary key A point having lower 
X coordinate value has higher priority of being scanned early 

• Amom^st the points with equal X coordinates, those representing right 
end points have higher priority of being scanned early (This ensures 
that at a particular position of the scanline the deletion of edges has 
higher priority than the insertion of edges) 

• Sort the group of points that have the same X coordinates and all 
points representing the left end points in increasing order of Y 
coordinates 

Step 3 Sort the elements of VfJ array according to the following strategies 
• Sort the points in increasing order of Y coordinates 
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• Amongst the points ^\■lth equal Y coordinates the top end points have 

higher prionh than the bottom end points 

• Sort the group of points with equal Y coordinates, and all representing 
the bottom end points in increasing order of X coordinates 

Step 4 Maintain an arra\ V_sc[ ] that represents a vertical scanline The 
elements oj V_sc lepresent the points that are currently on the scanline 

Step 0 Run the scan line ^_sc[ ] from left to right From H[ J array select a 
group oJ points that have equal X coordinates 
Label 5 1 for (each point e group) 

Begin 

if (point = left end point) 

Be^in 

Inset t that point in the scanline V_sc[ ] at proper position such 
that the elements of V_sc[ ] are always in the y-sorted order 
End if 

else replace the corresponding left end point by the current right 
end point 

End for 

step 6 Between every pair of end points P1,P2 e V_sc[ J such that 
(PI repieseiits a valid width constraint, join a constraint of weight 

equal to the minuniim width of the corresponding layer Delete all right end 
points from the scanline V_sc[ ] 

if (all groups are not selected) select next group and goto Label 5 1 

Step? Maintain an array H_sc[ ] that represents a horizontal scanline The 
elements ofH_sc represent the points that are currently on the scanline 
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Step 8 Run the scan line H_sc[ ] from bottom to the top From HI J array 
select a group of points that have equal Y coordinates 
Label 8 1 for (each point e group) 

Begin 

if (point = bottom end point) 

Begin 

Inseri that point in the scanlme H_sc[ ] at proper position, such that 
the elements of H_sc[ ] are always in the x_sorted order 
End if 

else leplacL the corresponding bottom end point by the current top end 
point 
End for 

Step 9 Between every pair of end points P1,P2 e H_sc[ ] such that 
(PI P2 }},„„) lepreseiits a \alid width constraint join a constraint of weight 
equal to the minimum width of the corresponding layer Delete all top end 
points from the scanlme H_sc[ } 

if (all groups are not selected) select next group and goto Label 8 1 


2 3 Spacing Constraint Generation 


The generation ot spacing constraints is a very time consummg process m 
mask la\out compaction based on constramt graph approach In worst case there are O 
(Is ) constraints where N is the no of pnmitives in the lavout One wa\ to limit the 
numbei ot constraints is not to generate redundant constraints Shadow propagation is 
one tiaditional approach to generate the spacing constraints where X MSibihty and Y 
visibility are checked to determine whether two edges of layout entities are to be 
constrained bv the spacing constraints Two primitives are said to be Y visible if a 
hoiizontal line extended Irom one pnmitive in +Y direction meets the other without 
passing through tn\ extra pnmitive Synonymously Y visibility can be defined 
Concept ot Y visibilit\ is shown m Figure 10 Here the pnmitive 5 is Y visible from 
the piimitixe -1 while primitive C and D are not 



figure 10 Illustrating the concept of visibility Here J?isJYvmMefrom-4 while C and 
D are not 
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The conbtraint generation method used should ideally generate an irredundant 
set ot constraints since the cost of sohmg the constraint graph is proportional to the 
number of edges ot the constraint graph In practice we must trade mcreased 
constraint generation time to generate the smaller constramt set for smaller soKing 
time 



F igurc 1 1 Illustiating Shadowing using X Visibility 


The shadow propagation ilgorithm examines ihe positions of the cells to 
dcteimine what constraints are redundant Refemng to Figure 11 we can imagine a 
legion that contains the rectangles that must be constrained against primitive A in 
horizontal direction as tailing under a shadow cast from 4 If the shadow falls on a 
lectangle that rectangle must be constramed against A Hence we must generate 
constraints from 4 \o B C and D We need not generate a constramt firom^ to £* as £" 
does not fall in the shadow of A However E is indirectly being constramed agamst \ 
through the entit\ C as E falls under shadow of C which is in the shadow of^ 

The abo\e mentioned method is good enough for generating a irredundant set 
of constraints but the only problem with this approach is that it produces an under 
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estimated set of constraints The diagonally interactive entities cannot be constrained 
b) this approach hence this cannot prevent two initially non o\erlapping primitives 
trom overlapping each other in new technologies Figure 12 is an illustration of this 
tact 



4 



(a) 


(b) 


[ iguie 12 IllubUating the shortcomings of shadowing using X-visibihty (a) initial 
iwontiguration beloie compaction (b) final configuration after compaction 


In figuie 12(a) the initial layout with all its rectangles spaced from each other 
lie shown The arrow represents the vertical and horizontal spacing constraints 
between the vaiious rectangles as generated by conventional shadow propagation 
ilgoritlim The lectangles A and D are imtially non overlapping However if the 
layout IS compacted as there is not any constraint between A and D these two 
primitives cannot be prevented from being overlapped as shown m Figure 12(b) To 
ensure the spacing relations m diagonal direction we modify the defimtion of shadow 
and redefine it bv casting a shadow at an angle of 4^ as shown m Figure 13(a) 
Shadowing m di igonal direction ensures that the primitives m diagonal direction, e g 
A and D are constrained and when the layout is compacted the topology can be 
preserved The basic idea of not generatmg redundant constramts is as follows 
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Figuic 13 (a) lllusuating concept of diagonal shadowing Here A and D are 

constrained (b) I he la\out attei compaction Since .4 and D are constrained, the} will 
not overlap 


It a primiin c C lies in the shadow of another pnmitne B which is also in the 
shadow ot the pnmiine 1 then a spacing constraint from 41 to C is redundant even if 
C lemains in the shadow of .4 However this statement holds good, if the lavout is 
maae up of boxes onlv In practice the basic layout pnmitiv es are not only boxes, but 
also the polvgons The difficulty with the polvgons as basic layout entities is 
illustrated in Figuie 14 Heie B and C are in shadow of^ and C is in the shadow of 
B however the constraint (A-^C) is not a redundant one. and actual redundant 
constraint is (B-^C). Fuithci in the cases of rectangles, if there is a spacing 
constraint from piimitive/4 to primitive C. there cannot be a spacing constraint from 
C to 4. which IS not valid m actual cases with polygons and boxes as the basic layout 
primitives In Figure 14, both the pnmitives ^ and C are related in both wavs 
Keeping all these aspects in view we propose and validate an algorithm, that does not 
uenerate any redundant spacing constraint. Instead of shadow being defined by the 
primitives, we define shadovv trom every right edge of the primitive while 
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constraining in horizontal direction and from the every top edge of the pnmiti\e while 
constrainina in \ emcal direction Thus we need the followins definitions 

Shadow. An edge is said to be in the shadow of another edge if it is 
diagonally visible from ej 

Front- An edge is said to be in front of another edge g/, m honzontal direcnon if 
the following conditions are satisfied 

• Both ei and e 2 belong to the same pnmitive 

• Edge tTj IS X visible from the edge g; 

Lejt Edge- An edge e, of pohgon P is said to be a left edge, if the left side of e, is 
\ acant and right side of edge e, is filled to an extent, at least the minimum resolution 
ol the layout Relerring to Figure 1 5(a), the edge e, ePcan be called a left edge of i* 

Right Edge- An edge e, of poh gon P is said to be a right edge if the right side of the 
edge IS vacant and the left side is filled to an extent at least the mimmum resolution of 
the la\ out In Figuie 1 5(a) the edges ej„ e^^, Sjs are the right edges of the polygon P 

The concept of shadow and front is illustrated in Figure 15 In Figure 15 (a), 
the edges e„ ej„ ej 2 and belong to the same pnmitive and and are X 

visible trom e, Thus e^j, ej 2 , and ej^ are said to be m front of e^ According to the 
above stiategies, in Figure 15(b) edge e,? is vsx front of edge e,/, ^j3 front 

ol edge ej,. and e,^ is m front of e,j The edges e^j is diagonally visible from the edge 
e,. So edge ej, can be said to be in the shadow of e^ Similarly e^ is m the shadow of 
e,-! and eji Synonymously top edge and bottom edge can be defined to consider 
shadow in the vertical direction 
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2.3.1 Checking for Left and Right Edge 


The concepts of left and nght edges of an entir\' are defined earlier To 
determine whether an edge e, e P, is a left edge of a pnmiti\e. we ha\e to predict a 
point to the left or nght of the edge and check whether that point belongs to the 
polygon Two subcases are descnbed below Here we predict a pomt P(x,y) to the 
nght of the edge e,, and use point enclosure algonthm to check whether this pomt is 
inside of the polygon P 




Figure 1 6 Check for left and right edge of a pnmitive (a) case of a vertical edge, 
(b) case of an inclined edge 


CASE A: e, is a vertical edge : 

Let P(x,y) be a test point sufficiently close to the edge e, towards its right as shown in 
Figure 16 (a) The coordinates of the point be 

X^X, + 5 
Y= (T/ + F/;/2. 

where Jis a small real number, smaller than the minimum resolution of the layout 
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If the test point P(x,y) remains inside the poHgon. then the edge e, must be a left edge 
of the polygon, otherwise it is a right edge. 

CASE B. e, is an inclined edae: 

S' 

This case is shown m Figure 16(b) In this case the coordmates of the test 
point P(x,y) will be 

X=(X! +X^)/2 + S, 

Y=(Y! + Y:)/2. 

If the test point is inside the poh gon. then e, is a left edge, otherwise it is a right edge 


2.3 2 Checking for Front 

The concept of front has been defined earlier Here a method to check whether 
an edge Cj is va front of another edge gj is described Before the edges are checked for 
being front, the minimum width type constraints are generated All edges e, that 
remains m front of e, are constrained by e, Thus checking for front can be achieved 
by simply checking for adjacency between the corresponding nodes in the current 
constraint graph .An outline of the procedure Front( ) is given below 

Procedure Front(e„ej) 

Begin: 

e, eP, and ej ePj; 

V, = the node in constraint graph corresponding to edge ej 
Vj = the node in constraint graph corresponding to edge ej 
if (P, ^Pj) return (FALSE) ; 
if (!Adjacent(v^Vjj) return(FALSE) ; 
return(TRUE); 

End; 
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2.3.3 Checking for Diagonal Visibility 


The range of the coordinates m XY plane for which an edge is said to be 
diagonally visible from another edge e, is decided by the onentanon of the edge m the 
layout plane Here we will consider only the construction of horrzontal constraints and 
hence, we will restnct ourselves discussing the diagonal visibility m honzontal and 
diagonal direction only hence we limit the orientation of an edae to be either vertical 
or inclined onl\ For different orientation of edges e^ and the limits of diagonal 
visibility IS descnbed below In all cases we presume e, to be a right edge and e, to be 
a left edge 

CASE 1. Both e, and e, are vertical lines: 

Figure 17(a) shows a general case when ej is diagonally visible from e, 
Figures 17(b) and (c) show the two extreme limits for diagonal visibility We define 
two variables and F/^^ that define the two extreme limits for the pair (e„ e^) of 
edges Y/„gi, and Y/g^ are determined to be 



Figure 17 Illustratmg the diagonal visibility Case of both vertical edges 
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Y,., = Y,‘.(X,-X,) 

where F," is the F coordinate of the top end point of edge e, and F/ is the F 
coordinate of bottom end pomt of the edge Similarly F/ and Y/ are defined If 
^ ^ ^/iigh or Yj < Yio^ then Bj does not remams in the shadow of e, or Cj is not 
diagonally visible from 



(a) (b) 


Figure 1 8 Illustrating diagonal visibility. Case of g, being an mclined edge 

CASE2 : e, is vertical and Cj is inclined line: 

The two sub cases for this type onentation are shown m Figures 18(a) and (b) 
Let X, be the X coordmate of the lower end point of Cj and X 2 be the X coordinate of 
upper end point of For both the sub cases the extreme values of the F coordinates 
are also shown in Figures 18(a) and (b). Yhigh ^low take the values* 
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^high Y, + (Xj- Xj) 

Yio^ =Y/-(X2-X,). 


Yj > Yf„g/j ox Yj < Yig^, then ej is not diagonally visible from 
the shadow of 


else Cj remams m 



Figure 19 Illustrating diagonal visibility Case of e, being an incimed edge. 

CASE 3: is inclined line, ej may be vertical or inclined line: 

This case is shown in Figure 19 In this case, for simp licit}', we restnct the 
definition of diagonal visibility, and define the shadow considenng only the 
honzontal visibility. We obtam the Imuts as 

Yhigh ~ Y^ and 

Y =Y^ 

■* low I 
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2.3.4 Checking the Shadow 


Using the above concepts of left edge and right edge, and the diagonal 
visibility vve can formally define the shadow as follows 

Shadow An edge e^ is said to be in the shadow of another edse e, if 

• e,€P, and e, is a right edge of P, 

• e^e Pj, and e, is a left edge of Pj 

• gj IS diagonally visible firom e, 

A formal algorithm to check the shadow is outlined below 

Procedure H_Shadow (e,^ Cj) 

Begin: 

If (ej lies to the left of e,) return(FALSE); 

If (Cj IS a right edge) return (FALSE); 

If (e, IS a left edge) return (FALSE); 

= Upper limit of shadow region in Y direction for the pair of edges <e„ ej>. 

Yiow - Lower limit of shadow region in Y direction for the pair of edges <e„ ej>. 

If (Y bottom of ej > Yh,gh) return (FALSE); 

If (Y top of ej < YiJ return (FALSE); 
return (TRUE); 

End; 


2.3 5 Generating the Spacing Constraints 


For generating the spacing constraints, we propose a method analogous to 
depth first traversal m a directed acyclic graph [13],[14].[15]. The edges of 
primitives of the layout resemble the nodes m the DAG, and we traverse along the 
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pnmitive edges one by one depending on the front and the shadow cntena as the 
DAG IS tra\ersed from one node to the other along its arcs The algonthm starts from 
a source vertex (S) that represents the left boundaiy^ of the la\out and visits ever} 
edge before stopping The algonthm is cited below 

Procedure Generate_H_Space_Constramts(Gfj) 

Begin: 

Stepl .’Create two arrays of edge listing called Leftf } array and Rigltt[ ] array. 

(a) Left array contains all the edges that represent a left edge of a pnmitne 
and the right side boundary of the layout (T),ue. the elements of Left[ ] array 
can be the destination vertex of a spacing constraint. 

(b) Similarly the Right[ / array contains all the edges that represent a right 
edge of a primitive and the left side boundary (S) of the layout, ue., the 
elements ofRightf ] can be the source vertex of a spacing constraint. 

(c) Sort both Leftf ] and Rightf ] arrays according to the following strategies: 

. Sort the elements in increasing order of X coordinates. 

. For a set of elements with equal X coordinates, sort in increasing order 

t^f ^bottorrr 

. The first element of Leftf ] array is the edge corresponding to the left 
boundary (S). 

. The first element of Rightf ] array is the edge corresponding to the 
right boundary. 

Step 2: Create an array Vmled[ I that maintains the visa status of an edge. 
Inttiahze all of its elements to FALSE. VisitedfOJ to Visitedtno__ofJeft_edges - 
U are attributed for keeping the visit status of edges in Leftf f array, and 

Visitedfno_ofjefi_edgesj to risited[no_ofjeft_edges + no_of_right_edges/ 

keep the visit status of the edges in the Rightf ] array. 

Step 3: root = S; 
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for (each edge e, e Left[ J) 

Begin: 

if (I Visited (e)) HJTraverse (root, e) ; 

End for; 

End; 

Here H_Traverse( ) recursne routine for traversing the edge listings, and other 
procedures such as Min_space( ) and Descendant_Visited( ) are outlined below 

Procedure H_Traverse(p, q) 

Begin: 

If (q IS a left edge) 

Begin: 

Visited(q) = TRUE; 

Min_space = Minimum required spacing between edges p and q; 

If (p and qfall in interactive layers) join constraint from p to q of value 
equals to Mmjspace; 

for (each edge ej in front of node q) H_Traverse(q, e) ; 

End if; 

else If (q IS a right edge) 

Begin: 

Visited (q) = TRUE; 

for (each edge e* ^ Left[l 1 is in shadow of q) 

Begin: 

If(IVisited(eJ) H_Traverse(q, e^ ; 

else if(Visited(eO and e* h not visited by any descendant of q) 
Join spacing constraint from q to e^ of weight equals to the 
minimum spacing between q and 
End for; 

End if; 
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End; 


Procedure MmJSpace(p,q) 

Begin • 

Ciirr_entry = first entry of design rule table of spacing rules; 
wliile(all entries are not searches) 

Begin: 

If ((Curr_entry — flayer 1 =p -flayer) and 
(Curr_entry —flayer! = q —> layer)) 
return (Currjentry Spacing_value) ; 

Curr_entry = next entry of the design rule table. 

End while; 
return(-l) 

(-1 corresponds to non interactive layer) 

End; 

Procedure Descendant_Visited (p,q) 

(Returns true if node p is visited from any of the descendant of node q) 
Begin : 

If (Adjacent (p,q)) return(TRUE); 
for (each node n, adjacent to p) 
if (Descendant_Visited(n„q)) return(TRUE); 
return(FALSE); 

End; 
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2.4 Generation of Overlap/Extension Constraints 


0\erlap constraints abstract the overlap '’extension lavout rules These rules are 
imposed to preserve the interconnection among the lav out elements in different 
lavers/levels Vanous overlap 'extension rules have been descnbed in an earlier 
section on design rule descnption In this section we will describe the vanous 
algorithms that have been developed as a part of this thesis work to generate these 
constraints A basic scanline approach is emploved to detect the overlap among the 
lav out pnmitives Special care is taken to generate constraints corresponding to the 
contacts and vias The following is an illustration of the scanlme approach for 
generating overlap/extension constraint 

For generating overlap extension constraints in honzontal direction, i.e. to 
construct a honzontal scanline is run in the v ertical direction on the layout plane 
Tlie flattened layout plane information is stored in a common data base, e g . an array 
containing the vertical edge listings This arrav is sorted according to a certain 
strategy such that when the scanline is run. at anv instant (or position of the scanline) 
some selected pairs of edges will generate overlap constraint Whenev er a bottom end 
point is encountered, that edge is inserted to the scanline and when a top end point is 
encountered, the corresponding edge is deleted from the scanline. The number of 
distinct steps or positions that a scanline can take is equal to the number of distmct Y- 
values of the end points At a particular position of the scanline an edge can be 
inserted or deleted and the deletion takes higher pnonty than the insertion. From 
among a set of edges falling on the scanlme. a valid overlap is checked according to 
the following strategies For a particular position of the scanlme, let and ej be the 
two edges cut by the scanline The vanous conditions under which an overlap 
constraint between e^ and ej is generated are 


45 



• If the edges are not pre\iously constrauned If these are pre\ioush 
constrained, then it is either b\ minimum width tv'pe constraint which 
implies both e, and Bj are the edges of the same primiti\ e (polygon or box) 
or It IS b\ a spacing constraint which implies these edges should be 
spaced 

• If there exists a valid design rule which specifies the required overlap 
benveen the lav out elements For example in Figure 20. the overlap 
between e,i and Cj; is valid while the overlap between e,j and ej 2 is not valid 
as there is no such design mle ( CM to CO 

• If £?, eP, and Cj P, and Pj are different pnmitiv es and P^n Pj # e g . 
for an overlap constraint to be generated between two edges, the 
corresponding pnmitives must ov erlap. i e one pnmitive mav be 
completel} embedded within the other 



Figure 20 Illustrating existence of valid overlap rule 


So at a particular position of the scanline we select all the edge pairs satisfying 
the above criteria and constrain them The most vntal problem that can arise with this 
approach is the case of contacts and vias as described below; 
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In Figure 21(a), a contact CC and a metal Ia\er CM is shown. If the scanlme 
approach is used, in this case the following constraints will be generated. 

in Gh and 

®m5 ^'“c 4 and e^i ^ejjj 3 in Gy 

After the corresponding L P problem is solved, the contact position in the 
resultant la\out will be as shown m Figure 21 (b) At this position it will \iolate the 
design rule between and e^j, and between and e^g as there is no such constraint 
between this pair of edges ongmally generated 

One way to eliminate this problem is to generate constraint between e„i and 
2 nd between e„g and However with these constraints the wire width will be 
unnecessanh wider Instead if we generate constraints 2 nd (e„^—>€^i) we 

prevent the contact position being shifted into the active wire, hence can set the wire 
width to Its minimum width. However these constraints cannot be captured by the 
conventional scanlme approach since the edge leaves the scanlme before enters 
and e^i lea\es the scanlme before enters So we treat the contacts and Mas as a 
special case and generate the constraints pertinent to the contacts and vias pnor to 
running the scanlme 
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Figure 21 Illustrating shortcoming of scanlme approach in generating constramts for 
contacts and vias (a) Imtial contact position, (b) Finally contact is shifted mto the 
wire, violatmg design rule. 
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Figure 22 Illustrating constraint generation for contacts and Mas Blow out box is 
shown in dotted line Edges which are constrained are shown b> arrow 

For adding the constraints pertinent to the contacts and Mas. we take one 
contact or Ma at a time We then determine the layout pnmimes which are actually 
connected b\ this contact For example, we have to search for two n_diffusion/ 
p_diffimon/ poly/ metall polygons which encloses that contact. Similarly for a Ma we 
ha\e to search a metall polygon a metall polygon which enclose that Ma. Once such 
a tnplet, e g . CM-CC-CP, is detected, we blow out a rectangular region of some 
arbitrary large length and width surrounding the contacu'via. The dimensions of the 
blow out region depend upon the minimum required overlap between contact and the 
characteristic layer This is illustrated in Figure 22 in which a metal wire and d. poly 
wire are connected through a contact The corresponding blow out box is shown in 
dotted line All the edges that fall inside or intersect with the blow out box are taken 
into consideration for generating constramts We can ensure constraints between 

(em6-> ^p 2 > and (ern 2 ->^p:i) by this method. 
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2.4.1 The Algorithm 


Procedure GeneratejOverlapjConstraint () 

Begin • 

Step 1: Generate Constraints for contacts and vias; 

Step 2: Create an array H[ J whose elements are the end points of all 
horizontal and inclined edges in the layout Create an array V[ J whose 

elements are the end points of all the vertical and inclined edges in the 
layout 

Step 3 Sort the elements of H[ J array according to the following strategies: 

• Primary key: increasing order ofX coordinates. 

. Amongst a group of points with equal X coordinates those which 
represent right end points have higher priority. 

. Sort the group of left end points with equal X coordinates in 
increasing order of Y coordinates. 

Step 4: Sort the elements of V[ ] array according to equivalent strategies 

Step 5: Create an array V_sc[ J that represents a vertical scanline, whose 
elements correspond to the points that are currently on the scanline. 

Step 6: Run the scanline from left to right 

Label 6.1: From H[ ] array, select a group of points that have equal X 
coordinates. 

for(each point belongs to Group) 

Begin: 

If(point = left end point) Insert that point in the scanline V_sc[ ] at 
proper position, such that the elements of V_sc[ J are in the Y sorted 
order. 
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If(point right endpoint) Delete the corresponding left end point from 
the scanline V_sc[ ]. 

End for; 

Step 7: for (every point p^ eVjscfJ) 

Begin: 

for(every point Pj e H_scf ]) 

Begin: 

e, — the edge of the layout corresponding to p^; 

Cj = the edge of the layout corresponding to p^ 

Vj = vertex in in Gy corresponding to e„ 

Vj = vertex in Gy corresponding to 
P, IS the primitive containing e^; 

Pj IS the primitive containing e^; 

If (! Adjacent(v„v^) 

If (Valid_Overlap(v,yj)) 

If (Min_Overlap(v^Vj) ^1) 

Join constraint between v, and Vj of weight equals to 
Min_Overlap(v„Vj); 

End for; 

End for; 

If (all groups are not selected) goto Label 6.1: 

Step 8: Create an array H scf } representing the horizontal scanline. 

Step 9: Run the horizontal scanline from bottom to the top. 

Label 9.1: From V[ ] array, select a group of points with equal Y 
coordinates. 

for(each point belongs to Group) 

Begin: 
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If (point ~ bottom endpoint) Insert that point in the scanline H_sc[ J at 
proper position, such that the elements of H_sc[ / are in the X sorted 
order. 

If (point — top end point) Delete the corresponding left end point from 
the scanline H_sc[ ]; 

End for; 

Step 10: for(every point p, e H_sc[ ]) 

Begin: 

for(every point Pj e H scf }) 

Begin: 

e, = the edge of the layout corresponding to p,; 

Cj = the edge of the layout corresponding to p^; 

V, = vertex in in corresponding to e,; 

Vj = vertex m Gh corresponding to e^; 

P, IS the primitive containing e,; 

Pj is the primitive containing e^; 

Ij (I Adjacent(v„v)) 

If (Valid_Overlap(Vi,Vj)) 

If (Min_Overlap(v„Vj) ^1) 

Join constraint between vi and vj of weight equals to 
MinjOverlap(v,yj) ; 

End for; 

Endfor; 

If (all groups are not selected) goto Label 9.1: 

End; 


The minirQUiii overlsp requirements between, various lasers are stored in a linked 
list To find out the minimum overlap between two nodes, the list is searched to find a 
match of the layers Following is an algorithm for the function Min_Overlap( ) - 

central l<brari 
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Procedure Min_Overlap(v„Vj) 

Begin * 

layer 1 =v^—^ layerjtype; 
layer 2 = Vj ^^y^t‘_type; 

If (both V, and v, are left edges) lay er2 = -layer2; 
else If (both v, and Vj are right edge) layer 1 = -layer2: 

(In other two condmons, ue.l: v, « left edge and v, is a right edge then 
layerl and layer! will not be changed. 2: v, a right edge and v, left edge 
implies a spacing constraint and will not arise in this case.) 

Search along the design rule list representing the overlap rules for a match 
of layerl and layer2. If such a match found, return the corresponding 
overlap value of that entry. 

If such a match not found return (-1). 

End; 


n 


.5 Transistor Resizing 


Resizing of all the transistor present in the onginal la\out can be achieved by 
two ditterent ways The first and easiest approach is to scale all the transistor 
uniformly setting new length at some user specified value and maintaining the same 
W/L ratio as the previous technology The other approach is to scale the transistors 
according to user specified values Although the first approach is simpler, the second 
approach provides the user more flexibility in mampulatmg transistor sizes for the 
new technology We provide both the options The various types of transistor 
configurations that our tool can support are as shown in Figure 23 
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(a) ( b) ( c) 

Figure 23 Illustrating vanous transistor configurations (a) Simple transistor (b) 
Rectilinear Transistor (c) Transistor with 45° inclined feature 


The transistor resizing information are stored in a table whose elements are 

. A point P(x,y) in the layout plane, specify mg the position of the transistor, 
to identify a particular transistor This point corresponds to the S_W comer 
of the transistor 

. W/L ratio for the transistor in old technology 

. Required W/L ratio of the transistor in new technology 

For constant W/L scaling the only information required from the user side is 
the minimum channel length of the transistor The channel length and width wall be 
set to a value determined by the W/L ratio of the transistor in old technology For 
resizing the transistor according to new W/L, the first task is to detect all the edge 
pairs which have to be constrained Figure 24 illustrates the vanous cases that may 
arise In Figure 24(a), a simple transistor configuration is shown and the required 
constraint and its weights are shown In Figures 24(b) and (c) the required constraints 
for other configurations are shown We propose the following formal algonthm to 

resize the transistors 
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Procedure Resize_Transistor( ) 

Begin: 

for (each transistor in the layout) 

Begin: 

Stepl: Find the diffusion polygon Pdiff and the poly polygon Ppoly 
defining the transistor using the point enclosure algorithm . 

Stepl: From the polygons Pdiff and Ppoly determine the relevant edges 
that define the actual transistor area. 

StepS: Constrain the appropriate edge pairs according to the new W/L 
values, as shown in Figure 24. 

End for; 

End; 



Figure 24 Showing the vanous required constraints for vanous transistor 
configurations 
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2.6 Other User Specified Layout Constraints 


Although all the general types of constraints are discussed, there arc still some 
user specified topological constraints, known as user constraints (UC), which ha\e to 
be imposed to provide some flevibility to the user in retargetting the cells The \anous 
user constraints that our tool can support are 

• Fixing power (V^p and bus width to a user specified \alue 

• Fixing the connector positions in the la\ out 

• Fixing the cell height 

2.6 1 Fixing the Power Bus Width 

For any layout the power (Vpp and buses carry' the maximum current, 
hence, these have to be kept wider than the other metal interconnection wires 
Furthermore, in standard cell based design, the power buses bems common to all the 
cells located in the same row have to be of a fixed width 

For fixing the Vpp and bus widths, these buses have to be identified at 
early stages The connectors connecting the Vpp and buses are identified m input 
processing stage This case can be considered as a special case of generating 
minimum width type constraints, hence can be taken care of at the stage of generating 
minimum width type constraints In this case a pnmitive given as an input to the 
procedure Add_Min_Width_Constraint is first checked for its membership to the 
metall layer and then is checked for its being a power bus from the information of the 
power bus connectors If it is a power bus, the weight by which the vanous nodes 
corresponding to that entity are constrained is the mmimum width of a power bus 
rather than the minimum width of a metal layer intercormecL 


55 



2 6.2 Fixing the Cell Height 


For the standard cell based design, all the cells m the same row, have to be 
of the same height, hence it is essential that the cell height be fixed to a user specified 
value Cell height can be fixed by explicitly joimng a constraint between the two 
extreme nodes, i e , source (Sy) and sink (Ty), of Gy of weight equal to the required 
cell height 


2 6 3 Fixing Connector Positions 

Keeping cell abutment in view, there must be a common platform for 
positioning the I/O pins of a cells For this reason the connectors of a cell ha\e to be 
fixed at some position along the cell boundary A standard rule is to fix the connectors 
at a distance of (integer + 1/2) of the routing pitch, from the ongin of the cell (S_W 
comer of Abutment box) along the cell boundary All the cells designed according to 
the above strategy have their connector positions fixed along the cell boundar>'. hence 
can be abutted easily The connector positions can be fixed by identifying the 
corresponding mask polygons of met(ill/metcil2 layer, the connector is connected to, 
and joining explicitly constraints from source to the respective polygon edge of the 
value equals to the distance of the connector from the source 
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Chapter 3 


Solving the Constraint Graph 


3.1 Introduction 


In Chapter 2 we discussed the process of generating constraint graph, in hich 
the layout problem is formulated as a linear programming (L P) problem The L P 
problem consists of a set of constraints and an objecti\e function The aim is to solve 
this problem in order to optimize the objective function, keepma in view that all the 
constraints are satisfied Furthermore, in our case each equation in the L P problem 
has at most two variables Hence the problem in algebraic form is translated into a 
graph theoretic form The following is an illustration 
Let the set of constraints be as follows 

Xb-Xa > 8 1 

Xc-Xa > 12 2 

Xc-Xb >10 3 

Xd-Xb>10 4 

Xd-Xc> 5 5 

Xe-Xc>10 6 

Xe-Xd> 2 7 

These set of constraints are represented by a directed graph as shown in Figure 25 
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F igure 25 Illustrating graph representation of a set of linear equations 


The task at hand is to select a group of constraints that leads to a soluuon to 
the problem Among the set of constraints some may be redundant, e a . constraints 2, 
4 and 6 m Figure 25 are redundant constraints Our task is to detect and remote all 
such redundant constraints present in the constraint graph 

It is possible that there may be several groups of constraints and that each 
group satisfies all the design rules There anses the problem of choosing a group that 
minimizes certain parameters specific to the latout, which is called optimization 
Each group of constraints leads to a different basic feasible solution. We follow the 
Simplex method for solving the linear programming problem, but the problem is 
actually solved in graph domain rather m conventional tabular manner 

The key idea of simplex method in solving the linear programming problem is 
to select a basic feasible solution (bfs). This imtial bfs may not correspond to an 
optimal solution From this bfs, successively, one or more constraints are removed 
and equal number of constramts are added without violating any design rule, and m a 
manner to reduce the objective function value In this chapter we discuss the 
algorithms and various relevant issues related to findmg a basic feasible solution to 
the problem, formulation the objective function and finally we will discuss algonthms 
and issues related to the optimization of a layout 
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3.2 Basic Feasible Solution 


Formally linear programming problem can be defined as follows 

Instance : A real m xn constraint matnx A, a real m-\ector b. and a real 
«-vector c 

Configuration: All real n-vectors 

Solution. All /i-vectors x that satisfy the linear constraints A x>Q 
Minimize C x where C= cost vector 

The sets of feasible solutions to the constraints when represented in R“ plane, 
define a region (may be bounded or unbounded) called a PoKtope E\erv point inside 
the polytope is a solution to the problem The basic task of optimization is to find one 
among those points that minimize the cost function Although eveiy' point inside the 
poly tope represents a feasible solution, the following lemma restncts the search 
region of the polytope and minimizes the number of points to be checked for 
obtaining the minimum cost function 

Lemma : 

Let P be the polytope defined by the set of inequality constraints and c(x) = C^x 
be the linear cost function Then the minimum cost on P with respect to cfv) is 
attained m a comer of P 

So instead of checking for the objective function value at any arbitrarv' point 
inside the convex polytope, the check is restncted only to the comer points of the 
poiytope When represented in graph domain each of the comer points of the poKiope 
corresponds ro a spanning tree Each basic vanable corresponds to a node m the 
constraint graph (also in the spanmng tree) Each inequality' (facet of the poKiope) 
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corresponds to a constraint m the constraint graph Conventional approach of sohing 
linear programming utilizes a tabular format to store all the information But the 
la>out problem consists of a large number of vanables for which the tabular format 
for stonng and solving is quite space expensive Comentional approach for sohing 
L P problem by simplex method is descnbed in [16] Here, we descnbe onl> the 
algorithms and methods to solve the L P problem using the graph theoretic approach 

Conceptually our basic objective is to obtain a layout of mimmum overall 
area To achieve this the circuit elements and \ertical interconnections are moved 
horizontally to the leftmost achievable position Similarh the horizontal 
interconnections are moved to the bottommost possible location Presenth we 
generate vertical and horizontal constraint graphs taking all the lower bound type 
constraints into consideration, i e all the constraints are of the form 

X.-X, >Wt. 

where Wt is a positive real number specifying the required mi mm um separation 
between the two nodes Thus our constraint graph G is in fact a single source single 
sink, edge weighted directed acyclic graph 

We define a solution set of G as the set { V/x), 0 <i <n } such that it satisfies 
all the inequality constraints implied by the edges of G Since Vq represents the source 
(SJ, eg, the left boundary of the layout, its position is normalized to zero A basic 
feasible solution set { V^(x) / of G is defined as the minimum solution set, if, for any 
other solution set {Vi(x)} we have V^(x) < V^(x) All the nodes are placed at as 
minimal distance as possible from the source node The value of VJx), e g , distance 
of sink vertex is the value of mimmum possible width of the compacted layout. As the 
primary goal is to minimize overall area of the layout, all the nodes that fall in the 
longest path from the source to the sink remain invariant to the optimization process. 
The only nodes that can be moved durmg optimization phase are those which do not 
fall on the longest path from the source to the sink and which provide some slacks for 
movement 
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3.2.1 Removing Redundancy from the Constraint Graphs 

The constraint graph so far constructed ma> have a number of redundant 
constraints Some of these redundant constraints can be removed before determimng a 
bts to the problem The removal of redundant constraints helps to speed up the 
subsequent optimization steps To remove all redundant constraints requires the 
following graph reduction For each constraint .Y, - X, > the longest path length 
from K, to F, is calculated If the longest path length is larger than then the 
constraint is redundant and should be removed. The task of graph reduction is 
equivalent to the construction of the transitive closure of the graph, which is though 
polynomial in time is still too expensive O(N^) Instead vve use the simpler task of 
removing those arcs which satisfy triangle inequality Let us consider three vertices 
K are constrained by weights and as shown in the Figure 

26 We will remove the arc (V, V) if the tnangle inequality^ + ^kj holds 

good Following is a formal algonthm to trim some of the redundant arcs from the 
constraint graph G 



(a) 



(b) 


Figure 26 Illustrating removal of redundant constraint Here < ff^i^+fVkr 
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Procedure Arc trimmins 
Begin 

for (each arc e,j = (Vi,Vj) belongs to G) do 
Begin 

for (each arc = (Vi,Vk) belongs to G) do 
Begin 

if (arc Cj^j = (Vk, Vj) belongs to G exists and 
Wij <Wik + mj) 

Remove arc e,j (\%Vj) from the graph G; 

End: 

End; 

End; 

The number ot possible arcs going out from a node F, or is on an average, a 
constant, independent of the size of N. since most of the arcs are between neighbonng 
nodes The average time compiexitv' of the arc tnmming algonthm is some constant 
times l^l The worst case complexity can be as bad as that of fmQing the transitive 
closure 

3.2.2 Searching a Basic Feasible Solution 

After a selected set of constraints is removed from the constraint graph, we 
solve the longest path problem m order to find a basic feasible solution. The following 
IS an algorithm for determining the bfs to the problem 

Procedure Searchjbfs 
Begin 

1. From the graph G select a spanning tree T. 

2. Compute the distances of all nodes of T. 
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3. Check whether ,he tree Tthus selected leads to a feasible solutwn. 

4. If the tree is a feasible tree return, 
else 

Make an elementary tree transformation. 

goto 2. 

End 


The abo\e procedure starts from an arbitrar> spaniung tree selected from the 
constraint graph The tree thus selected may not satisfy ail the design rules The edses 
ot G which are to be added to the spanning tree in order to make it feasible are 
detected and equal no ot edges are deleted from the spanmng tree bv a process called 
tree transiormation I he tree transtormations are repeated until there are no infeasible 
arc present in the resultant tree This procedure is applicable to both the honzontal and 
venical constraint graphs to determine a bfs The steps defined m the above algonihm 
are explained and the algonthms for these steps are desenbed in the following section. 

3 2 2 1 Selecting Spanning Tree from the Graph 

The spanning tree problem is basically an edge selection problem More 
precisely, given an edge weighted graph G = (V,E), the problem is to select a subset 
of edges E’cE such that E' induces a tree T = (V,E) that spans every vertex of the 
set V The algorithm to select such a spanmng tree is desenbed as follows 

Let G be a directed graph, free of selfloops (in fact the constramt graphs for 
the layout problems are self loop free) Let G consists of n \ertices and e edges. The 
vertices are labeled 1,2,... n and the graph is represented as a linked list. For selectmg 
a spanning tree, the edges of the graph are selected one by one and mserted into the 
tree if it does not produce any loop in the resultmg tree Each vertex is assigned to a 
set Thus at start there are n different sets At any stage each set represents a partial 
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spanning subtree constructed so far For a connected graph C, at the end. all die sets 

are merged into a single set ivhich contains all the vertices and a selected no of edges 
defining the spanning tree 

Lemma : 

In the constraint graph (7, there is at least one directed path from the source 
node to ever>' other node There is at least one directed path from every node to the 
sink node 

In the process ot selecting the spanmng tree, the root of the tree corresponds to 
the source node and one of the leaf node corresponds to the sink node We are 

required to select a subset of edges such that there is one and only one 

directed path existing from the source node to every other node This implies that 
e\ er\’ node m the spanning tree must have one and only one m degree, the in degree of 
the source node being zero 

Define two arrays D[ ] and S[ Jos 

For a node u belonging to the spanning tree, D[u] = TRUE if the in_degree of u is 
not equal to zero Tins implies the node u remains at the destination of some edge in 
the spanning tree Similarly for a node u belonging to the spanmng tree, S[ul = TRUE 
if the out_degree of node u is not equal to zero That means the node u is the Source 
of some edge belonging to the tree If at any stage, an edge (u,v) is inserted to the 
spanmng tree formed so far, D[v] becomes TRUE and it will prevent any other edge 
(u\y) to enter to the spanmng tree Thus all the edges for which D[v] is TRUE, are 
reiected from the spanning tree and these are marked as co_tree edges If D[v] is 
FALSE, an edge (u',v) can be inserted to the spanmng tree. But the exact method of 
insertion depends upon the value of DfuJ, SfvJ, and SfuJ Figure 27 is an illustration 
An algonthm to select a spanning tree from the graph is given below 
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(d)Oll 


(e) 100 


(f)101 





Figure 27. Trees arranged in order of D[u], S[v], and S[u]. Dotted edge (u,v) represents 
edge to be inserted to the trees Solid edges represent the edges of the trees formed so far 
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Procedure SeIect_Spanning_ Tree(G) 

Begin 

1. Create a tree node and assign its various parameters corresponding 
to the source node (SJ of graph G 

2. for (all node u) Initialize : 

D[u} FALSE; 

S[uf = FALSE, 

3. for (each edge (u,v) e G) do 
Begin: 

if(DM = TRUE) Begin ■ 

Mark edge e (u,v) as cojree edge ; 

Continue ; 

End if; 

Check Dluj,SM, andSfuJ; 

Case 000 : 

Neither ii nor v is included in any of the subtrees formed so far. Create two 
tree nodes u and v. Assign various parameters to u and v according to the 
Loirespondim’ node in graph G. Join an edge from u to v and mark v as the 
son oj u. 

Case 001 : 

Node It already exists in one subtree while node v is not in any of the 
subtrees. Create a tree node v, make suitable mapping between node v and 
corresponding graph node. Link u and v appropriately. 

Case 010 : 

Node u does not exists in any of the subtrees, while node v exists in one of 
the subtrees. Create a tree node u and make suitable mapping between u and 
corresponding graph node. Join an edge betsveen u and v such that v will be 
recognized as a son of u 
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Case Oil . 

Both the nodes u and v are present in separate subtrees Tu and Tv. Link u 
and V appropriately such that the Pvo subtrees Tu and Tv will be merged into 
a single tree. 

Case 100: 

In some subtree, there is an edge from some node u' to node u and node v 
does not exist. Create a tree node v. Join an edge between u and v to make v 
the son of u. 

Case 101 

Node u IS in one of the already constructed tree while node v does not exist. 
Create a tree node v and link u and v appropriately. 

Case 110 : 

Nodes u and v exist in two separate subtrees Tu and Tv. Join an edge from u 
to v so that V will be son of u and the two subtrees will be merged into a 
single tree. 


Case 111 : 

Nodes u and v exist in two separate subtrees Tu and Tv. Link u and v 
appropriately to merge the two subtrees into a single tree. 

End for ; 

End ; 
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3 2 2 2 Computing Distance of Every Node of the Spanning Tree 

After selecting a spanning tree firom the constraint graph, the next task is to 
compute the location of vanous nodes Smce the mformauon about the spanning tree 
IS stored in binary tree form, computing distance can be achie\ed bv a prs-ordsr 
tra\ersal of the binary tree The source node fS’J is fixed at the location '0' The 
location of successive nodes are determined by adding the edge weight to the distance 
of their father node For an example, let be the position of node edge eu^(u,v) 
belongs to tree T and = weight of edge e^,(u,v). Then the position of node v is 
gi\ en by Xy = X„ + W, The following is a recursive routine that tra\ erses the tree in 
pre-order and determines the position of all nodes 

Procedure Compute_dist(T) 

Begin: 

if(TI=NULL) 

Begin : 

if (T—> father — NULL) dist(T) = 0; 
else 

Dist(T) = Dist(T -> father) + edge weight; 

Compute_distance(T —>■ son) ; 

Compute_distance(T brother) ; 

End if; 

End; 

3 2 2 3 Determining Feasibility of Tree 

Tlie sparming tree selected so far may not lead to feasible solution, that 
satisfies all the design rules, i e , some of the constraints ongmalh present m the 
constraint graph may not be satisfied The following lemma suggests the condition 
under which the selected spanning tree represents a basic feasible solution. 
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L mma Every edge (u v )0 T defines one fundamental circmt/loop m the tree T 
D tme length of a fundamental loop defined by an arc (u v)eT as 


L ^ }Vt(v,vj) X Dir(ViV^ for all ^(v belong ng to the fandame t 1 loop 

\ her V]) is an arc in the fundament 1 loop defined bv the ed e ) 

Dir(ii v) = +1 if the edge (VjVj s in same direction as the direction of arc fr v) 

Ise Dir(u v)"l 

The cot responding spanning tree T leads to a fe sible solution if for all 
lundamcntal loops th length of the fu damental loop is non positi e In a formal 
wav we sav a tree T is feasible if for all cotree edges (uvje T L ^ < 0 

The following is an explanation of the above lemma 

Let us consider the spanning tree as sho vn in F gu e '’8 The dotted arc 
icprcsents the edge of the cotree and the solid a cs represent the ed es of th tr e 
Ed e weights are spe tfied along with the edges Consider an edge (u \) of the cotree 

The fundamental loop defined bv the edge (uv)\s u-->v 

The direction of edges a shown in Figure '’8(b) The length of the fundamental loop 
= 10 10 + 10 5 8 10 = 7 which implies that the onginal tree is not a feasible 

one The explanation for this mfeasibility is as folio vs 

The node v is not placed at its minimal possible position it is placed at 23 
units awav from node but if ve follow the path v -+ m -+ v we will obtain 
the position of node v as 30 units from node v which is the min mal (left most) 
position of node v If we restrict ourselves to the ongmal spanning tree it will lolate 
the design rule between the node u and node v These two nodes should be separated 
by a minimum distance of 10 umts while these are actually separated by onU units 
only To make the solution feasible we have to delete the edge (v v) from the ongmal 
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spannin tree and add edge (u v) to the tree This is called an element ry tree 
tr nsformation in which one spanninfi tree is extracted from anoth r on by emo mg 
one edge from the initial spanning tree and then adding another ed e to it Here e 
will discuss the algontlim to check, whether a gi en spannin tree leads to feasible 
solution If the tree T leads to an infeasible solution we vill det ct the exact d 
(ii\)eT which should be added to the spanning tree to make tlie soluuon feasible 



(a) 



Figure 28 Illustration of a tree leading to a feasible solution 
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Procedure Check Feas biIty(T) 

B III 

for (e ery d e (u v) eT) do 
B in 

t t(u y) = weight of arc(u vj 
L ^ - wt(u v) + Dist(u) Dist(v) 
if(L ^’>0) 

Be in 

Infeasible arc = (ii y) 
r ttirn(FALSE) 

Lnd if 
End for 
ret urn (TRUE) 

End 


There are E N +1 edges present in cotree Hence the time complexitv of the 
al orithm is 0(C) 

3-24 Tree Transformation 

In the previous stage we detected the edge (u v) which must be added to the 
spanning tree in order to make the solution feasible Simultaneously one of the edge 
\hich is originally present in the spanning tree must go out of the tree Anv edge ejj of 
ih fundamental loop that belongs to the T and with direction +l mav be a suitable 
candidate tor removal Without loss of generality we choose the edge (u vje T to be 
removed Following is the algorithm for tree transformation 

Procedure Tree_transformation(T u v) 

Be'^in 
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1 Detect the edg (a v) s T 

2 Mark d<fe(xi v) as cojtree ed e 

3 Mark td<fe (u v) as tree edge 

4 \Iak suitable modification in the on inal tree T to remove d e 
(u v) and to add the ed c (ii v) 

End 

3 3 Optimization 
3 j 1 Need of Optimization 


The bnsic feasible solution obtained i pre\ lous st p shnnks e% erv elem nt to 
a coiner (origin) of the layout plane and produces a layout of minimum overall area 
The overall area of the initial layout does not depend on all the pnmm es/elements 
onh a few control the layout area We call those elements as the ntical elements 
specific to a ^iven lavout instance However the noncntical elements have plenu ot 
freedom to position themselves in the original lavout within som 1 m ted ran e 
without affecting the overall layout area One ot the important considerations n 
mi rating the layouts is to preserve the uniform distnbution of la out elements 
throughout the entire layout plane The layout obtained from the solution of the 
longest path problem suffers from the following demerits 

\1I the non critical elements are packed towards one om r of the la out 
which introduces higher degree of non uniformity to the layout 

The resultant layout has a longer total wire length Hence it contnbutes to 
more parasitic resistance which degrades the circuit performance 
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Figures 29 and 0 illustrate the importance of wire length m n mization In 
Fi u e ^9(a) Uie blocks A B !md C are interconnected by vires and die total \ i e 
len th amounts to 1 units This configuration is obtained afte sol in^ the Ion est 
path problem and without wire length optimization Figure 29(b) is the correspondina 
la out obtained after uirelength optimization Here net wire I ngth = 9 units 



(a) (b) 


Fi ure 9 Illustrating wire length minimization 

In Figure 0(a) one end of the wire \V is pulled to the left most ach e\able 
position while its other end is connected to an element on the cntical path Hence the 
net wire len th is unnecessarily higher In Figure 0(1) th same layout is sho^vn 
after optimizing the wire length Here the net wire length is reduced and hence results 
m an improved circuit performance without compromising the net layout area In 
Figures 31 and 3'^ we present more cases where circuit performance can be improved 
by optimizing the layout Consider a contact location where two wires belonging to 
different lavers are interconnected together as shown m Figure 31 Here wirel is m de 
up of metall and wire2 is made up of polysillicon Resistivity of polv is orders of 
magnitude higher than that of metal Hence wires on polv la>er significantly 
increases si nal delay In Figure 31(a) the contact is placed at its mmimal possible 
position However the contact has some freedom m positive X direction In Figure 
31(b) the same layout is shown with the same net wire length but with the contact 
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shitted towards right This results in a decrease in vire length on poly layer while the 
vire length on the metal J layer increases by an equal amount B> this proc ss the 
o\ rail length of the wire remains the same as its previous value but the effective 
r sistance is reduced by significant amount vhich ensures bett r circuit pe formance 



r igure 30 (a) A layout compacted without wire length minimization 
(b) The same layout compacted with wire length mm mization 

In standard cell based design the I/O sp cifications hence the connector 
positions are fixed along the cell boundary Figure 32 shows a cell for which the 
connector position is fixed Consider the case where a inetal2 w e associated with a 
connector making contact with a metall wire (Figure 32(a)) In Figure 32(a) the 
lav out IS obtained after placing all the elements at their minimal possible positions 
vhich lesults m extending the metal2 wire and shrinking the metall wire If we 
assume metal2 layer is more costly than metall layer then the layout involves mo e 
cost F ure 3'^(b) shows an alternative arrangement for the same functionalitv but 
vith shorter wire length for metal2 layer This arrangement is a result of wirelength 
optimization and involves lesser cost 


74 






r%uic31 Moving the contact to reduce o erall parasitic resistan e 
(.a) Layout vithout vire length optimization involves hi her resistance 
(b) Layout after wire length optimization Length of wre n poly la r s 
d ci eased expandin wire on m tal layer Overall len th remains same but 
IfeUive r sisiance reduces significantly 


Y X 

■I 1 I 


« 4 





(a) (b) 

figure 32 Illustrating vire length optimization X represents connector position 
(a) Layout before optimization involves more wire length on metal2 (more costly) 
laNcr Hence involves more cost(b) Layout after vire length optimization Wire 
length in m fal2 layer is reduced by increasing wire length on metail layer Overall 
resistance is reduced 
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3 3 2 Formulating the Cost Function 


riie chaiacteristic parameter specific to a layer type which is of our nterest is 
the sheet resistance of the laver The sheet resistance of anous la ers are e^en n 
T ibl 4 Chapter 4 Various vays of formuiat ng the cost function are descnbed in 
[17] and [18] In [18] two different st at is ha\e been conside ed to model the cost 
finction One of them imoKes a relation vith the perimeter of the mask pol\ on to 
model the cost function vhile the other takes into account the area of the mask 
pol aon to model the cost function In our approach we associate ach ed e \ith a 
cost parameter which ve define a given below 

Definition 

For an arc (uv) associated with a particular laver the cost of arc C is 
d tilled bv the rate of change of resistance of the associated mask la er vith respect 
to the length of the arc If the arc is not associated with a single la er or not at all 
associated with a layer its cost is set to zero 

L t for an arc (u v) R = Resistance of the corresponding mask considering 
the direction (u->v) as the length of the sheet and the other dimension as the \ idth of 
the sheet Let/? be the sheet resistance of the corresponding mask laver 
L - Length of arc vj 
W = Width of the sheet (refer to Figure 33) 

CostC =rf/? /dL 

R =R X Number of squares on the sheet defined by arefw v) 

= R x(L /WJ 
Thus CostC =R /fy 
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As the distance between u and v changes while determinin the teasible 
solution len th and width of the sheet should refer to the dimensions of the most 
recently available sheet So we determine the cost parameters after etting a feasible 
solution to the problem The cost function to be optimized can be formulated as 



1 iguic 33 This polygon represents a particular mask layer Cost of arc (ii ) will be 
decided by the incremental resistance of the shaded region 

Z = X C x(X XJ + £ C x(Y YJ 

itveGfi iiveGi 

The two different terms are optimized separately in two passes one in horizontal and 
the other in vertical direction 

We have assigned a cost to each of the arcs which represents the incremental 
resistance with respect to the length of the arc Consider a node which is coiine ted 
through a number of such arcs Moving that node towards nght will reduce some arc 
length and increase the other The immediate conclusion is to move that node as much 
towards nght as possible if such movement results in a decrease in the total 
resistance We define the cost parameter for each node as folio \s 
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Defimtwn For a node v the cost of v (CJ is defined as the me ease m the 
r sistance per iiovement of the node v towards nght by unit distance 


So C IS the sum of the costs of arcs which connect the node to some other 

nodes i e C = 2^ (C C J for all nodes u connected to v So the cost fiinction for 
oplimuation is now linearized and can be iven as 

z- ZC xx 

Consider a component containing nodes V; v v If we mo this 
component is a whole by unit distance the net ncrease in cost will be qual to the 
mere isc in cost of indn idual nodes and can be given by 

C(T0= ECi 

Vj e 

3 3^ The Algorithm 

The optimization problem is solved using the simplex algorithm sp lalized tor 
Inndlin linear constraints modelled by graphs The input to the optimizat on problem 
IS the feasible solution set (feasible spanning tree) obtai d n the pre ous stage This 
feasible tree serves as the starting solution for simplex algonthm In th s section ve 
will describe the optimization procedure in detail Following is the algonthm 

Procedure OptimizeJLayout(T G) 

Be in 

1 Determine the Ion est path from the source node (S) to the sink node (T) 
Mark the nodes that fall on the longest path as critical nod s and the arcs on 
the longest path (critical path) as critical arcs 
for (each node v e G) 

Begin 
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if(v IS not a critical node) 

Be in 

Detenmne the cost of movm<^ that node by unit d stance towards ri bt 
End If 
End for 

for (each arc ej(v vJeT) 

Be in 

if(e j IS not a critical arc) 

Begin 

1 Partition the tree T into tivo subtrees T1 and Tl by s ici f! at v e T1 
and Tl contains the source node and v^eTl 

2 Coinpiii the cost of movi ig Tl sc whole by unit 1 stance vfiiclt is 
qiial to sum of costs of node v v&Tl 

3 if (cost < 0) 

B in 

Det rntine the arc of G whi h offers m ninti m slack for the novem nt 
of Tl Add that ed e to tl e original t ee aid remove ed^e (v v) from 
the tree Make suitabl modifications m the on mal graph G 
R CO npute the dista ices of va wits nodes 
End if 
End if 
End for 

End 


In the above algorithm we are considering one component at a time Each 
component is defined by one of the non cntical edge of the feasible tree If we remove 
a non critical edge e/v, vj) from the tree T the tree will be partitioned into two 
subtrees T, and Assume that the node v, belongs to the subtree T; which also 
contains the source node (S) and the node v; belongs to the other subtree (Ti) Note 
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that all the critical nodes will aK ays remain in the subt ee Tj The nodes belonain to 
T, will remiin undisturbed tor the current pass and the nodes belon in to Tj will 
exploit the possibility of movement to reduce the cost function 

Th cost of moving all the nodes in T/ b unit distance as a whole is 
consid red If this cost is nonpositive it impl es that it is beneficial to move the 
component The next step to be considered is the maxim 1 distance by which \ e can 
move this component Tj Tl e interaction of the set of nodes in Tj vith the rest of 
nodes vill decide this maximal distance The component 7/ mav have different 
bhcks vitli respect to the nodes m 7/ ho ever t is the minimum slack that n eds to 
b considcicd This is becai se ve cannot move a component by a distance more than 
Ih minimum ivailable slack The correspondn edge must b used to merge the two 
sublrc s T, and 7, into a single tree and the edee e ,(v vy is to be removed If the cost 
of moving the component as a whole is positive the ed e e / v) is not emoved from 
the initial tree but marked as \ isited and ill not be considered further When all the 
non critical edges are visited the algorithm terminates and the resultant tree obta ned 
corresponds to the optimal solution 

An example of optimization process IS sho vnmFi ures 4 35 and 6 Figure 4 

d scribes the initial constraint raph in hich the ed e veights and edge costs are 
sho vn The cost function IS gi enby Z-10(X6 \3)+10(X4 Xi) + 100(X7 V-#) 

Figure 35 represents one basic feasible solution to the problem The positions 
of various nodes as per this solution are sho vn The cost of this solution is calculated 
to be '>'^50 The longest path from source to sink is sho m in dashed line The nodes 
Vfl, V Vs V; and Vj are in the longest path hence are the cntical nodes The positions of 
these nodes are invariant to the optimization process The nodes v, Vj and have 
some freedom for movement towards right The cost of moving these nodes are C ; - 
IOC 3 ^ 10 C 4 = 90 mdC 6 = ^0 We can immediatel conclude that by moving 
nodes v, Vj and towards nght will reduce the cost Movmg will increase the cosh 
hence it will not be moved 
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fiLUie 36 Optimal solution to the problem The positions of vanous nodes are 
shown Cost of this solution = 650 


Although mov ng the node reduces the cost its movement s restn ted by 
the node But when the edge (v vj is considered the nodes ^ and ^ leman n 
som component and the effecti e cost of moving the component as a whole - 80 
and there is some slack for the movement The edge ( $Vs) provides the minimum 
slack Hence the edge (v^vg) is added and fvj is removed Figure 6 shows the tree 
corresponding to the optimal solution Tlie positions of vanous nodes are sho m The 
cost of optimal solution = 650 which shows significant reduction m the cost 
compared to the solution obtained by solving the longest path problem 
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Chapter 4 

Results Conclusions and Scope for Future Work 


4 1 Results 


In order to validate the applicability of the proposed algorithms all th 
mentioned algoritlims are impl mented using C progxammm Ian iiaae and run on 
UYI\ based HP 9000 workstation A number of expenments ha e been made on 
some of the existing cells of ARCUS T cfmology Limited \ ith a set of desien rules 
specified in [3] All the layouts are retareetted to a 0 5 micron process The anous 
mask layers and the detailed listing of all the design rule applicable for this process 
are given in [3] We took the sheet resistance (ohms/square) as the charactenstic 
parameters pertinent to a mask layer which ar provided to us by M/s ARCUS 
T chnolof^ Ltd Table 4 shows values of these parameters Table 4 1 sho vs the 
results obtained for different t st cases and Table 4 2 depicts the time taken by vanous 
steps for those test cases From these results it is evident that the constraint generation 
phase takes the maximum time and is as expected This is because of the extra 
computation necessary to generate as few irredundant constraints as possible vhich 
in turn speeds up the latter stages of layout recycling 
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Table 4 1 Performance analysis of TECHMIG for vanous test cases 


Pmamcteib of comparison 

Layout #1 

Layout #2 

layout TO 

No ol pnmitives 

9 

30 

39 

No of nodes in HCG 

3 

70 

126 

No of nodes in VCG 

23 

70 

126 

No of arcs m HCG 

38 

253 

488 

No of arcs m VCG 

40 

237 

563 

Input ceil size (p ’) 

11\6 

28v26 

13 8 8 

Output cell size 

8^\44 

12 9 \ 18 6 

12 1 X 19 6 

Cost lunction value before 

13698 

2393283 

4366063 

optimization 




Cost function value after 

3173 

2394110 

3760038 

optimization 





Table 4 Timm analysis of the various stages of TECH\tIG with respect to anous 
test cases shown m Table 4 1 


CPU time in second 

Input 

I/O 

CG 

Obtaining 

Wire Length 

Total 

Layout 

Processing 

Generation 

a b f s 

Optimization 

Time 

Cell # 1 

0 17 

0 03 

0 04 

0 01 

0 27 

Cell ' 

0 24 

0 98 

0 02 

0 03 

127 

Cell# 

oh 

8 01 

0 36 

0 13 

8 81 
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T ible 4 3 Sheet resistance of various mash Ia> ers 


Name of the layer 

Symbol of the layer 

Sheet resistance (Q / □) 

N diffusion 

LC\D 

3 0 

P diffusion 

LCPD 

2o 

N well 

Lcmv 

740 

P veil 

LCPW 

740 

Polvsillicon 

LCP 

28 

Contact 

LCC 

38 

Metal 1 

LC\t 

0 0n3 

Metal2 

LC\I2 

0 033 

Vn 

LCC2 

1 7 


4 2 Conclusions 


A core tool TECHMIG for basic t chnolo y migration has been developed 
bv using newer and extending e istmg aleonthms for ns vanous building bio ks Tlie 
constraint graph compaction method has been adopted for design ng the important 
compactor block in TECHMIG In this method the topological descnption of a given 
instance ot cell layout is first translated into a pair of constraint graphs called the 
horizontal constraint rapli and the vertical constraint g apli respectively These 
graphs serve as a model for formulating the compaction problem as a 1 near program 
Several algorithms based on computational geometry have been employed for the 
topological analysis of a given layout instance and for generating the two constraint 
graphs corresponding to irredundant constraints These have been discussed in details 
A graph theoretic version of the simplex method has been coded to solve the 
optimization problem with the minimum delay wire length as its objective function 
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This objective function is modelled as the weighted sum of the wire len ths 
TECHMIG has been tested on a fe \ cells from an existing cell of \f/s ARCUS 
r htioiooy Limited The retarcretted cells are obtained without any user intenention 
In each case tlie final cell is not onh smaller it also has the total ire leneth 
optimized for minimum wire delays 


4 3 Scope for Future Work 

We have developed a basic ompactor tool TECHMIG Further nha cements 
need to be carried out in TECHMIG to provide it with more features and flexib lit\ 
so that It can result in more compacted layouts vhi h ha e better timing 
pertorminccs The present vork can be extended to include a number ot ssues as 
given below 

1 Hierarch} Preservation Lar e industrial designs generally include hierarchical 
descriptions of cells To aid the design of larger cells the present vork can be 
extended to support hierar hical descriptions of the input cells The technoloej 
migration of such cells can again be done in an hierarchical manner 

2 Ai tomatic Jog Insertion/DeleUon In a pract cal layout there may exist fe v areas 
where a straight wire can be joa ed to reduce the overall area of the layout \lso there 
may exist some zigzag portions of a wire which can be straightened to improve the 
timing performances of the cell by reducing the net wire length without affecting the 
overall layout area Such considerations have not been included n the present work 

3 Cr tical Path Breaking In optimizing a layout TECHMIG maintains the cntical 
path as an invariant in the optimization phase The idea behind such a cntical path 
indentation is that the portions of the layout which corresponds to the cntical path are 
exactly minimum according to the current set of design mles and user constraints In 
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general a cell can be made smaller by mov ng th objects off the critical paths into 
non cntical areas thereby creating parallel cntical paths This feature needs to be 
included m TECHMIG 

4 Contact Optimization If certam portions of the target layout permits more number 
ot contacts than in the source lavout without increasing the overall area of the cell 
incorporatma such extra contacts can result in lower parasitic and improved current 
carrying capabilities Such a fa ility can be added to the existing tool for obtain ng a 
better lavout A roup of parallel contacts can be merged to ether to form a super 
contact ol a larger size In the final lavout this super contact can be sliced into an 
apptoprnte number of contacts of the desired size 

0 Phmn Plo ving cm be used interactively to rearrange the geometrv of a cell b 
compacting a sparse lavout or creating, space in a dense la\out to esult m a bette 
fabrication yield The user places a plo p (a horizontal or vertical line segment) into 
the lavout and gives a dir ction and distance the plow is to move Acco dine to the 
specifications the plow is then moved through the layout and the material b hind it is 
stretched (if necessary) and the material in front of it is compressed 

6 W tal Layer Translatw Present tre ds in the layout design phase of VLSI 
chips involves the use of larger number of metal layers for routing to reduce the 
0 erall cell size Thus it may become essential to migrate exist ng cell layouts to a 
new process with larger number of metal layers Unfortunately none of the 
commercially available tools for layout conversion or technology mieration 
incorporates this essential feature 
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